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DISCOVERY OF NARROW X-RAY ABSORPTION LINES FROM NGC 3783 
WITH THE CHANDRA HETGS 
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ABSTRACT 

We present the first grating-resolution X-ray spectra of the Seyfert 1 galaxy NGC 3783, obtained with the 
High Energy Transmission Grating Spectrometer on the Chandra X-ray Observatory. These spectra reveal many 
narrow absorption lines from the H-like and He-like ions of O, Ne, Mg, Si, S and Ar, as well as Fe XVII-Fe XXI 
L-shell lines. We have also identified several weak emission lines, mainly from O and Ne. The absorption lines 
are blueshifted by a mean velocity of as 440 ± 200 km s" 1 and are not resolved, indicating a velocity dispersion 
within the absorbing gas of a few hundred kms" 1 or less. We measure the lines' equivalent widths and compare 
them with the predictions of photoionization models. The best-fitting model has a microturbulence velocity of 
150 kms -1 and a hydrogen column density of 1.3 x 10 22 cm -2 . The measured blueshifts and inferred velocity 
dispersions of the X-ray absorption lines are consistent with those of the strongest UV absorption lines observed 
in this object. However, simple models that propose to strictly unify the X-ray and UV absorbers have difficulty 
explaining simultaneously the X-ray and UV absorption line strengths. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — galaxies: individual (NGC 3783) 
— X-rays: galaxies — techniques: spectroscopic 



1. INTRODUCTION 

X-ray spectra of Active Galactic Nuclei (AGNs) often show 
evidence for deep O VII (739 eV) and O VIII (871 eV) absorp- 
tion edges. The ionized gas component creating these edges 
is named the 'warm absorber' and is seen in most (>70%) 
Seyfert Is and some quasars. Theoretical models for the warm 
absorber suggest it should also be a source of X-ray emission 
(e.g., Netzer 1996) and absorption (e.g., Nicastro, Fiore, & Matt 
1999) lines. While the current X-ray data are consistent with 
the presence of such lines (e.g., George et al 1998a), it has been 
impossible to study them in any detail. 

The bright Seyfert 1 galaxy NGC 3783 (V « 13.5 mag, z = 
0.0097) has one of the strongest X-ray warm absorbers known 
with a column density of ionized gas of Nn w (1-2) x 10 22 cm" 2 . 
ASCA spectra also show warm absorber variability and excess 
flux around 600 eV which has been interpreted as emission 
lines from the warm absorber, particularly O VII 568 eV (e.g., 
George et al. 1998a, 1998b). UV spectra of NGC 3783 show 
intrinsic absorption features due to C IV, N V and H I. Two kine- 
matic absorption components are known with radial velocities 
of « -560 kms" 1 and w -1420 kms" 1 relative to the optical 
redshift from [O III]. The C IV absorption shows extreme vari- 
ability (e.g., Crenshaw et al. 1999 and references therein). It is 
not clear whether the UV and X-ray absorbing gas are the same, 
although some relation seems likely (e.g., Shields & Hamann 
1997). 

This letter describes the discovery of a large number of X-ray 
absorption lines in the Chandra spectrum of NGC 3783. In § 2 
we describe the observation and data analysis, and § 3 gives the 
preliminary interpretation of the results. 



2. OBSERVATION AND DATA ANALYSIS 

NGC 3783 was observed using the High Energy Transmis- 
sion Grating Spectrometer (HETGS; C. R. Canizares et al., in 
preparation) on the Chandra X-ray Observatory 3 on 2000 Jan- 
uary 21. The detector was the Advanced CCD Imaging Spec- 
trometer (ACIS; G. P. Garmire et al., in preparation). The total 
integration time was 56 ks, and the observation was continu- 
ous. The data were reduced with the standard pipeline by the 
Chandra X-ray Center (CXC) on 2000 January 26. We used the 
Chandra Interactive Analysis of Observations (CIAO) software 
to repeat the reduction process using an updated response ma- 
trix file. Our reduction produced somewhat less noisy spectra 
than those produced by the pipeline. 

The HETGS produces a zeroth order X-ray spectrum at the 
aim point on the CCD and higher order spectra which have 
much higher spectral resolution along the ACIS-S array. The 
higher order spectra are from two grating assemblies, the High 
Energy Grating (HEG) and Medium Energy Grating (MEG). 
Each grating assembly produces two spectra for each order. Or- 
der overlaps are discriminated by the intrinsic energy resolution 
of ACIS. The zeroth-order spectrum of NGC 3783 shows sub- 
stantial photon pile up, and we leave the detailed analysis of 
these data to a future paper. We briefly comment, however, 
on the presence of a weak (ss 10 photon) off-nuclear (w 34") 
source located at a 2 ooo = 1 l h 39 ra 01 ?1, <5 2 ooo = -37°44'52;'9, co- 
incident with one of the galaxy's spiral arms. This source may 
be a luminous (Lo.5-8 kev ~ 2 x 10 39 ergss" 1 ) X-ray binary in 
NGC 3783, although it could also be an unrelated background 
source. 

The first-order spectra from the HEG and MEG have S/N ra- 
tios of w 5 and w 2.5, respectively (at around 7 A). The higher 
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Rest Wavelength [A] 

FIG. 1 . — MEG first-order smoothed spectrum of NGC 3783. The two upper panels are binned by 0.005 A, and the bottom panel is binned by 0.02 A. Fluxes are 
corrected for Galactic absorption. Dots with error bars show the typical 2a statistical error at various wavelengths. The identified absorption features are marked. 
The few identified emission lines are marked with (r) for resonance, (i) for intercombination, and (f) for forbidden for the He-like ions, as well as La and L/3 for the 
H-like ions. 



order spectra had only a few (« 3) photons per bin and will 
not be presented here. In order to flux calibrate the spectra we 
used the CIAO software to produce Ancillary Response Files 
(ARFs). The flux calibration error is estimated to be <10% 
below 6 A and <20% above 6 A (CXC, private communica- 
tion). We also corrected the spectra for the Galactic absorption 
of Nh = 8.7 x 10 20 cm" 2 and the cosmological redshift. 

The two MEG first-order spectra agree well with each other. 
The dispersions of the MEG spectra are 0.005 A per bin, and the 
spectral resolution is w 0.023 A (approximately constant across 
the entire HETGS band). The two HEG spectra also agree well 
with each other and with the MEG spectra, though they have 
lower S/N. The HEG spectra have dispersions of 0.0025 A per 
bin, and the spectral resolution is about two times better than 
for the MEG spectra. In Fig. 1 we present the mean of the two 
MEG spectra. This mean spectrum was smoothed using a sim- 
ple 3-bin boxcar filter. Since above ^16 A there are <4 counts 
per 0.005 A bin, the spectrum presented in the bottom panel of 
Fig. 1 is binned to 0.02 A. 

The overall shape and flux of the X-ray spectrum are 
consistent with those found in previous X-ray observations 
of NGC 3783 (e.g., George et al. 1998b). The Chan- 
dra fluxes in several energy bands are also consistent with 
those from a simultaneous 4 ks ASCA spectrum to within 
« 10%. The total flux measured from the Chandra spectrum 
is (2.0±0.2)xl0- u ergs cm" 2 s" 1 in the 0.5-2 keV band and 
(4.2±0.9)xl0 -11 ergscnrV 1 in the 2-7 keV band. When 
plotting the modeled ASCA spectrum on top of the Chandra 
spectrum, we find general consistency in the overall shape as 
well as with the oxygen absorption edges. 

The most striking new features seen in the Chandra spec- 
trum of NGC 3783 are narrow X-ray absorption lines. In order 
to measure these spectral features we interpolated over all the 



absorption lines and fitted the resulting continuum with a cu- 
bic spline which serves as the continuum reference point for 
measuring the lines. In Table 1 we list the information about 
selected absorption and emission lines. In Fig. 1 we identify 
some absorption lines which are not listed in Table 1 due to 
brevity and lack of model predictions. There are also some ab- 
sorption and emission features which appear to be present but 
are not identified yet. We used multiple lines from the same ion 
to help with line identification. We checked for the existence 
of these lines in each of the MEG spectra (and also in the HEG 
spectra when possible) and measured them in the smoothed av- 
erage spectrum. 

The measured lines' Full Widths at Half Maximum 
(FWHMs) are consistent with the MEG resolution which is 
« 1300 kms" 1 at 5 A and « 500 kms" 1 at 15 A. Hence, the 
lines are not clearly resolved (see discussion in § 3.2). We were 
not able to put better constraints on the velocity field with the 
HEG due to the small S/N. 

3. DISCUSSION 

3.1. Absorption Lines, Column Density and Ionization 
Parameter 

The strongest observed absorption lines in the X-ray spec- 
trum of NGC 3783 are from the H-like and He-like ions of O, 
Ne, Mg, Si, S and Ar, as well as Fe XVII-Fe XXI L-shell lines 
(see Table 1 but note that some features represent blends, and 
the values given are the summed EWs). 

As explained above, the lines are not resolved, and the only 
firm conclusion is that, given a Gaussian profile, their FWHMs 
do not exceed a few hundred km s" 1 . Given this uncertainty, we 
have taken an approach that relies solely on the measured EWs 
of the lines. We assumed a certain microturbulence velocity in a 
single line-of-sight cloud and calculated a theoretical spectrum 
with the photoionization code ION99, the 1999 version of the 
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noN and Emission Lines From NGC 3783 





Predicted 


Velocity 


Predicted 


Measured 




A 


Shift 


EW 


EW 


Ion 


(A) 


(kms- 1 ) 


(mA) 


(mA) 


Absorption Lines 


o A. V 1 


A 77S 
4. / Z5 


^7 1 4-1 1 7 

—J / tit J 1 / 


n 7 


S H4-7 A 


o A. V 




^1^4-707 


S 7 
O.Z 


7 ^4-7 ^ 


Qi VTV 
A.1 V 


^717 
J. 11 1 


Q7H4-7S7 
— VZUztZfi / 


1 A 
1 .4 


1 n 1 4-1 ^ 


Si XIV 


6.181 


—534=1=242 


7 7 


17.3=b4.6 


Si XIII 


6.648 


-586±225 


18.4 


14.7=1=4.8 


MgXII 


7.106 


-379±21 1 


6.6 


6.2=1=2.2 


Me XI 


7.473 


-280=1=200 


8.6 


8.0=1=3.6 


MgXI 


7.851 


-382±191 


15.6 


9.7=1=3.1 


Mg XIP 


8.420 


-463±178 


>19.9 


25.6=t6.6 


Ms XI a 


9.170 


-261=1=163 


>25.9 


34.5±10.2 


NeX 


9.362 


-384±320 


6.2 


16.4±4.5 


Ne X b 


9.481 


— 474±158 


12.0 


18.6±4.9 


Ne X a 


9.709 


-216±154 


>15.7 


29.7±8.2 


Ne X a 


10.239 


-410±146 


>24.7 


19.8=1=5.5 


NeIX a 


11.001 


-518±190 


>20.1 


50.4±12.2 


Fe XVII 


11.250 


-186=1=133 


23.3 


15.7=1=4.3 


Fe XVIIP 


11.435 


-629±131 


>20.5 


21.0±14.2 


NeIX a 


11.547 


-701±129 


>49.3 


27.2±7.4 


NeX b 


12.132 


-667±123 


69.8 


49.2±19.8 


Fe XVII b 


12.263 


-391±122 


43.1 


51.9±12.8 


O VIII 


14.635 


-20±102 


34.3 


24.8±12.2 


VIII 


14.821 


-182=1=101 


38.9 


27.8±12.8 


O VIII 


15.176 


-79=1=197 


42.8 


35.8=1=11.2 


Fe XVII 


15.264 


-471=1=196 


42.7 


30.6±8.4 


Fe XVIII 


15.628 


-624±249 


18.2 


27.7±7.9 


O VHP 


16.007 


-346±121 


>52.3 


43.9±14.2 


VII 


18.628 


-515±322 


52.1 


41.9±14.8 


N VII 


19.120 


-510=1=109 


11.2 


18.8±6.8 


Emission Lines 


MgXII 


9.300 


306=1=225 


15 


11.0=1=6.8 


NeX 


12.132 


444±123 


36 


44=1=24 


NelX 


13.447 


64±111 


19 


99=1=88 


NelX 


13.553 


464±110 


13 


163=1=140 


NelX 


13.699 


26=1=175 


41 


245=1=215 


O VIII 


18.967 


415=1=126 


460 


335=1=310 


N VII 


20.912 


161=1=100 


28 


118=1=92 


O VII 


21.602 


54=1=69 


99 


197=1=169 


O VII 


21.807 


229=1=68 


52 


135=1=109 


O VII 


22.102 


154±67 


201 


631=1=583 



0.08 



0.06 



CO 
> 

'=3 

o- 

LU 

"D 
CD 



CD 



0.04 



0.02 




0.02 0.04 0.06 0.08 

Modeled Equivalent Width [A] 

FIG. 2. — Measured versus modeled absorption line EWs. Filled circles 
are model EWs that are a sum of all lines contributing to the blend. Right- 
pointing triangles are lower limits for model EWs (see § 2 for details). The 
agreement between the model and the measured points is noticeable (Spear- 
man rank-order correlation probability of as 10~ 6 ). A line with a slope of unity 
is drawn to guide the eye. 

TABLE 2 
Best Model Column Densities 



"Blended with Fe absorption lines which are not taken into account in the model. 
b Blended with Fe absorption lines which are taken into account in the model. 
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code ION (Netzer 1996). We then compared the results with the 
observations. The most important ingredients of the modelare 
the microturbulence velocity, the column density, the spectral 
energy distribution (SED), the ionization parameter, and the gas 
composition. 

The atomic data required for calculating line optical depths 
are taken from various sources. The oscillator strengths for the 
H-like and He-like transitions of all metals are well known (e.g., 
Porquet & Dubau 2000 and references therein). This is not the 
case for the various iron L-shell lines where there is no standard 
published data set. In this paper we rely on various publications 
(e.g., Mason et al. 1979; Bhatia & Mason 1980; Cornille et al. 
1992, 1994; Phillips et al. 1996; Saba et al. 1999) as well as on 
the National Institute of Standards and Technology (NIST) data 
set. These were combined with the Phillips et al. (1999) line 
lists, which are based on solar X-ray spectra and served as our 
major source for the wavelengths of the iron lines. Only 3-5 of 
the strongest iron lines, per ion, have been used in the analysis. 

Our models are based on those in George et al. (1998b) that 
gave satisfactory fits to the 1993 and 1996 ASCA observations 
of the source. They include a weak UV-bump AGN SED (the 
weak IR case in Netzer 1996) with To. 1-50 kev= 1-75, a single 
constant density cloud with Nh = 1-3 x 10 22 cm" 2 , solar com- 



position, and a 0.1-10 keV ionization parameter of U x = 0.13. 
The standard (E > 13.6 eV) ionization parameter for this con- 
tinuum is 19.3 and the mean temperature is 1 .9 x 10 5 K. We have 
assumed Doppler profiles with V 2 = V^ ermal + V 2 rbulence and cal- 
culated absorption EWs for 100 < Vmrbuience < 600 kms" 1 . The 
best agreement with observations is for Vturbuience=150 kms" 1 . 
The model EWs are tabulated in Table 1. For a few lines we 
list only the model's EW lower-limit since the measurement in- 
cludes blends for which we do not yet have a modeled EW. 
In most cases, the blends are expected to increase the modeled 
EW by w 10-20%. The comparison between the observed and 
calculated EWs is shown in Fig. 2, and the deduced column 
densities, for the ions producing the strongest lines, are given 
in Table 2. All of the strong lines predicted by the model are ob- 
served in the spectra. The microturbulence velocity of our best 
fitted model is consistent with the velocity dispersion of the UV 
absorption lines (161±22kms _1 ; Crenshaw et al. 1999). 

We have carried out basic tests to verify that the deduced 
line widths are in the right range. If the 'true' widths are much 
larger than assumed, the resonance lines in the H-like and He- 
like series will be optically thin. The observed EWs show that 
this is not the case. On the other hand, if the 'true' line widths 
are much smaller than 150 km s" 1 , several of these lines will be 
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saturated, giving very similar EWs to all. This is not observed 
either. While these are not strong conclusions given the S/N of 
the observations, it is encouraging to see the good agreement 
between the observed and calculated EWs, suggesting that the 
line optical depths are indeed in the intermediate range. Ac- 
cording to the model, the largest optical depth for a line with 
E> IkeVis r(NeX 12.132 A)=35. 

3.2. Dynamics of the Absorbing Gas 

All the identified absorption lines in the HETG spectrum of 
NGC3783 are blueshifted by a few hundred kms" 1 with re- 
spect to the systemic velocity. The unweighted mean blueshift 
is 440±200 km s" 1 . Kaastra et al. (2000) have recently ob- 
served the bright Seyfert 1 NGC5548 using the Low Energy 
Transmission Grating Spectrometer (LETGS) on Chandra. The 
low-energy X-ray lines reported in that paper indicate a similar 
level of ionization, and the absorption lines are blueshifted by 
similar amounts. Judging by these first two results, it seems 
that the warm X-ray material in Seyfert 1 galaxies is outflow- 
ing with typical velocities of several hundred km s" 1 . This has 
important implications for the dynamics of the absorbing gas as 
well as for the X-ray-UV connection. 

The acceleration of X-ray ionized gas by the strong central 
radiation source in AGN has not been studied in detail so far 
(for some discussion see Reynolds & Fabian 1995). The recent 
work by Chelouche & Netzer (2000) is the first detailed study 
of this kind. This study discusses the motion of discrete, X-ray 
ionized clouds in hydrostatic equilibrium, under a large range 
of conditions. It shows that acceleration to velocities of a few 
hundred to a few thousand kms" 1 is a natural consequence of 
the AGN environment. The terminal velocity is of the same 
order as the escape velocity in the region where the flow orig- 
inates. Using these results, and assuming that the conditions 
in the nucleus of NGC 3783 are similar to those considered by 
Chelouche & Netzer (2000), we conclude that the origin of the 
warm absorbing gas in this source is outside the Broad Line 
Region, where the escape velocity drops below w 1000 km s" 1 . 
Better line profile analysis will enable us to constrain this loca- 
tion more accurately. 

Another important comparison is with the UV absorption 
lines observed from this source. Several papers (e.g., Mathur, 
Elvis, & Wilkes 1995; Shields & Hamann 1997 and references 
therein) suggest that the UV and X-ray absorbers may originate 
in the same component. The Crenshaw et al. (1999) observa- 
tions reveal two UV absorption systems, the strongest of which 



has a velocity consistent with the mean observed velocity of the 
X-ray gas and a velocity dispersion consistent with our best fit- 
ted model. Our model, which depends on the assumed SED, 
predicts column densities for C IV and N V (see Table 2) which 
are 37 and 5 times smaller, respectively, than what has been cal- 
culated by Crenshaw et al. (1999) using the 1993-1995 data. 
Hence the UV absorption lines appear to be inconsistent with 
our X-ray absorber model and are likely coming from a separate 
component. 

3.3. Emission Lines and Covering Factor 

The HETG spectra of NGC 3783 revel several weak emis- 
sion lines (see Table 1) some of which show hints of P Cygni 
profiles. These have a marginally significant redshift with an 
unweighted mean of 230±170 kms" 1 relative to the systemic 
velocity measured from the [O III] line. Kaastra et al. (2000) 
found similar redshifts for the emission lines in the LETG spec- 
trum of NGC 5548. While the EWs of some of the lines are well 
within the MEG and HEG capability, the low S/N and bumpy 
continuum do not allow reliable EW measurements. The val- 
ues given in Table 1 are therefore highly uncertain. However, 
even these rough measurements are good enough to test the idea 
that the central source in NGC 3783 is surrounded by gas clouds 
with properties similar to those deduced for the line-of-sight ab- 
sorber. This is done by comparing the EWs calculated by ION 
with those observed, for various assumed covering factors. 

The calculated emission-line EWs, for an assumed geometry 
of full covering, are listed in Table 1 . Given these values, we 
find that the mean deduced covering factor is close to unity, in 
agreement with previous estimates (e.g., George et al. 1998b). 
This is the first direct constraint on the covering factor of a 
warm absorber in any AGN. Unfortunately, the uncertainty on 
the covering factor is large, due to the uncertain EWs. Further 
discussion of this issue, as well as detailed analysis of the con- 
tinuum, emission features, and absorption features, are deferred 
to a future paper. 
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